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THE  VIEW,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT 
ARE  THOSE  OF  THE  AUTHORS(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS 
AN  OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION,  POLICY,  OR  DECI¬ 
SION,  UNLESS  SO  DESIGNATED  BY  OTHER  DOCUMENTATION. 


Research  was  carried  out  mainly  in  areas  of  time  series  analysis  and  multivariate 
analysis.  The  most  important  results  in  the  first  area  apply  to  regression,  autoregres¬ 
sion,  moving  averages,  and  more  generally,  stationary  stochastic  processes.  In  the 
second  area  the  results  have  to  do  with  elliptically  contoured  distributions. 

A  simple  condition  has  been  derived  for  a  numerical- valued  Markov  chain  to  be 
first-order  autoregressive,  i'lie  stationary  probabilities  of  a  general  finite  chain  can 
be  efficiently  estimated  by  the  mean  of  the  corresponding  vector- valued  process. 

Statistical  inference  about  regression  models  and  autoregression  models  is  often 
based  on  the  estimators  of  the  coefficients  being  approximately  normally  distributed 
-  at  least  in  large  samples.  'I'he  applicability  of  inference  based  on  zisymptotic  theory 
h2is  been  subsU;ii,;.-iiy  c.xtcnded;  a  Lindeberg-type  condition  is  sufficient  to  justify 
this.  The  disturbance  vectors  need  not  have  a  common  covariance  matrix  nor  do 
they  need  to  be  independent;  they  can  be  martingale  increments.  The  condition  is 
implied  by  a  strongly  uniformly  integrable  condition. 

A  very  general  asymptotic  distribution  of  a  (finite)  set  of  autocorrelations  has 
been  obtained  biised  on  the  theory  for  autoregression  coefficients.  The  results  for 
autocorrelations  are  valid  foi-  linear  processes  with  martingale  difference  disturbances. 

Goodness-of-fit  tests  foe  time  series  models  based  on  standardized  spectral  dis¬ 
tributions  have  been  developed  for  testing  patterns  of  dependence  in  stationary  pro¬ 
cesses.  Tests  can  be  conducted  on  the  basis  of  the  difference  between  the  sample 
standardized  spectral  distribution  and  hypothesized  process  standardized  spectral 
distribution.  Considered  a.'  a  stochastic  process  with  frequency  as  the  “time”  pa¬ 
rameter,  this  difference  multiplied  by  the  square  root  of  the  sample  size  converges 
weakly  to  a  Gaussian  stoclu.stic  procession  on  [0,  tt]  with  expected  vedue  0.  Because 
the  process  depends  on  autocorrelations,  the  conditions  for  the  limit  results  are  very 
weak.  For  a  linear  process  (that  is,  Anderson  has  shown  the  limiting 

distribution  of  the  autocorrelations  for  the  disturbances  being  martingale  differences. 

The  Gaussian  stochastic  process  is  characterized  by  its  covariance  function.  Af¬ 
ter  a  suitable  monotonic  transformation  of  the  frequency  parameter,  the  covariance 
function  of  the  transformed  process  is  a  constant  times  min(«,  u)  —  uv-j-  q(u)q{v),  0  < 
u,v  <  1,  where  q{xi)  depends  on  the  hypothetical  spectral  distribution.  Note  that 
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min(«,u)  —  uv  is  the  covariance  function  of  the  Brownian  bridge. 

Tests  are  based  on  functionals  of  the  process,  such  as  the  Cramer-von  Mises 
statistic,  the  Kolmogorov- Smirnov  statistic,  and  the  Anderson-Darling  statistic.  The 
distribution  of  the  Cramer-von  Mises  or  Anderson-Darling  statistic  depends  on  the 
eigenvalues  of  an  integral  equation  involving  the  covariance  function  of  the  limiting 
process;  the  integral  equation  can  be  converted  to  a  second-order  differential  equa¬ 
tion.  A  formal  solution  to  this  problem  for  standardized  spectral  distributions  has 
been  found,  although  it  is  feasible  only  in  certain  Ccises.  The  distribution  of  the 
Kolmogorov-Smirnov  criterion  satisfies  a  partial  differential  equation;  a  numerical 
solution  is  provided  by  solving  recursively  an  approximating  difference  equation. 

This  asymptotic  theory  has  been  jujlified  under  very  general  conditions  because 
the  distributions  of  the  criteria  depend  on  the  asymptotic  distributions  of  autocorre¬ 
lations. 

The  principal  investigator  with  Raul  P.  Mentz  studied  iterative  methods  for  es¬ 
timating  the  parameters  of  a  Gaussian  moving  average  process  of  order  1.  These 
procedures  yield  the  exact  maximum  likelihood  estimates.  Efficient  calculation  is 
based  on  the  calculation  of  quadratic  forms  of  relevant  matrices  in  the  time  domain 
and  the  traces  of  related  matrices  in  the  frequency  domeiin.  In  iterative  procedures  for 
exact  maximum  likelihood  estimation  in  first-order  Gaussian  moving  average  model 
“Woodbury’s  formula”  gives  an  efficient  method  for  evaluating  quadratic  forms  in 
the  inverse  of  a  covariance  matrix.  The  computation  of  the  tr2w:e  of  the  inverse  of 
a  covariance  matrix  can  be  expressed  as  an  infinite  series  in  the  autocorrelation.  A 
very  accurate  approximation  has  been  obtained. 

The  other  topic  of  research  is  inference  in  elliptically  contoured  distributions.  It 
was  shown  that  a  scale-invariant  vector-valued  function  has  the  saune  asymptotic 
distribution  as  in  the  normal  case  except  for  a  factor  depending  on  the  kurtosis  of  the 
distribution;  this  factor  can  be  consistently  estimated.  A  scale-invariant  scalar- valued 
function  with  value  and  first  derivatives  0  at  the  true  parameter  value  hais  the  same 
asymptotic  distribution  as  in  the  normal  case  except  for  the  kurtosis  factor.  In  most 
cases  this  asymptotic  distribution  is  x^- 
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